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ABSTRACT 

We investigate the thermodynamic stability of neutral black objects with (at least) 
two angular momenta. We use the quasilocal formalism to compute the grand canonical 
potential and show that the doubly spinning black ring is thermodynamically unstable. We 
consider the thermodynamic instabilities of ultra-spinning black objects and point out a 
subtle relation between the microcanonical and grand canonical ensembles. We also find 
the location of the black string/membrane phases of doubly spinning black objects. 



1 E-mail: dumitru@aei.mpg.de 
2 E-mail: maria.rodriguezQaei .mpg.de 
E-mail: theisen@aei.mpg.de 



Contents 

1 Introduction 



i 



2 Stress tensor and conserved charges 4 

2.1 Quasilocal formalism and conserved charges 4 

2.2 Doubly spinning solutions 5 

2.2.1 Black hole 5 

2.2.2 Black ring 6 

2.2.3 Black membrane [8 

3 Thermodynamic instability of the black ring 

4 Instabilities from thermodynamics |^ 

12 
15 



4.1 Ultra-spinning black holes . . 

4.2 Membrane phase of black rings 



5 Discussion 17 



A Temperature and angular velocities [20 



B Conditions for thermodynamic stability |21 



1 Introduction 

The physics of event horizons in higher-dimensional General Relativity (GR) is an inter- 
esting area of research not just for its intrinsic relevance to string theory. An investigation 
of black hole solutions in higher dimensions is also important because it has revealed new 
features: a richer rotation dynamics and the existence of regular extended black hole solu- 
tions. 

It is clear by now that some of the remarkable properties of four-dimensional black holes 
do not hold in general. A notorious example of particular importance concerns their horizon 
topology. In four dimensions, the spherical S 2 topology is the only allowed horizon topol- 
ogy for asymptotically flat (AF) stationary black holes. A related result is the 'uniqueness 
theorem', which states that a stationary AF vacuum black hole in four dimensions is char- 
acterized by its mass and angular momentum and has no other independent characteristic 
(hair) . 

The spectrum of stationary black objects is far richer in dimensions bigger than four 
(see [T] for a concise review). The restrictions on the topology of AF black holes require 
that spatial sections of the event horizon must be positive Yamabe type [2] and if spinning 
(stationary) they have to be axisymmetric [3] . The most obvious indication is the existence 
of a neutral AF solution describing a spinning black ring in five dimensions [HE]. As it 
was shown in jl], there are three solutions with the same asymptotic conserved charges (the 
same mass and angular momentum). On top of the well known Myers-Perry (MP) black 
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hole [6] with an S 3 horizon there are two different black rings with an S 1 x S 2 horizon 
(S 1 x S D ~ 3 in D dimensions). Therefore, unlike in four dimensions, the black objects in 
higher dimensions are not completely determined by a few conserved asymptotic charges0 
But the space of solutions of Einstein's equations in D dimensions also includes extended 
black holes. The black p-branes [8], dubbed black strings when p = 1 or otherwise mem- 
branes, are transverse asymptotically flat (only AF in D — p directions) , evade the no-hair 
theorems, and exhibit horizon topologies S D ~ 2 ~ P x MP. Interestingly enough, due to the 
richer rotational dynamics, in certain regimes in higher dimensions the thermodynamical 
properties of the compact black holes resemble those of the extended black objects. In this 
paper we aim to make progress towards a better understanding on these properties of black 
objects. 

We first study in detail the thermodynamic instabilities of doubly spinning black objects. 
We identify the ultra-spinning regimes (in parameter space) from their thermodynamic 
quantities. Our motivation stems from the observation made in [9] that, in dimensions 
greater than four, the thermodynamics of certain fast spinning MP black holes show a 
qualitative change in behaviour. That is, there is a transition towards a black membrane- 
like behaviour. This is due to the fact that, as one increases the angular momentum, the 
temperature of these ultra-spinning black holes reaches a minimum and then starts to grow 
as expected for the black membrane @ 

Black rings also exhibit a change in the thermodynamic behaviour that resembles the 
ultra-spinning regime of black holes when ultra-spinning along the S 1 direction. In the 
'thin ring approximation', when the radius of S l is much larger than the radius of S D ~ 3 , 
the singly spinning black ring can be approximated by a boosted black string. An inter- 
esting question we would like to address in this paper is how the ultra-spinning regime of 
the D = 5 black ring is affected by adding the second angular momentum, along the S 2 . 
In other words, whether neutral doubly spinning black ring is thermodynamically stable in 
the grand canonical ensemble and whether there is any connection with its ultra-spinning 
regime. 

Since there is no known background subtraction method for these black objects, we use 
a slightly modified version of the quasilocal formalism of Brown and York [TJj to address 
these questions. Supplemented with counterterms |l3j . the quasilocal formalism becomes 
a very powerful tool to study the thermodynamics of black objects that are AF. Recently, 
several concrete five- dimensional examples were discussed in detail in [7J [H] and in [15] for 
the unbalanced black ring. 

The relevant thermodynamic potential in the grand canonical ensemble is the Gibbs 
potential that is the Euclidean action divided by the periodicity of the Euclidean time [To] . 
The Euclidean method was applied to the black ring thermodynamics in [17J. Since the 
black ring does not have a real non-singular Euclidean section, the 'quasi- Euclidean' method 
[T8] was adopted to analyze the black ring thermodynamics. In this approach, the horizon 



*In [7] it was proposed that the necessary information to distinguish between black objects with different 
horizon topologies is encoded in the sub-leading terms in the boundary stress tensor. 

2 The numerical evidence of |10j suggests that the onset of the ultra-spinning regime for singly spinning 
MP black holes corresponds to a zero mode associated to the Gregory-Laflamme instability [11] . 
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is described by the 'bolt' in a complexified Euclidean geometry rather than a real oneo It 
was also pointed out in [T7] that the neutral black ring with one angular momentum is 
unstable to angular fluctuations — a more detailed analysis can be found in [7J [TH1 I2D] • 

By employing this method to compute the Gibbs potential, we find the response func- 
tions directly in the grand canonical ensemble. We observe that the second angular mo- 
mentum changes the situation in the sense that, unlike the black ring with one angular 
momentum, the doubly spinning ring is stable against perturbations in the angular velocity 
in some specific region of the parameter spacejj However, a careful analysis of all response 
functions that characterize the system reveals that the doubly spinning black ring is ther- 
modynamically unstable in the grand canonical ensemble. That is, there is no region in the 
parameter space in which all response functions are positive definite. 

On general physical grounds it is expected that the microcanonical ensemble of asymp- 
totically flat black holes is dominated by diffuse radiation states rather than black hole 
states [21]. In other words, it is favorable for the black hole to decay away (the heat ca- 
pacity is negative — see Appendix B for a detailed discussion on local thermodynamic 
conditions) and so pure thermal radiation is a local equilibrium state. Indeed, it was shown 
in [22] that, for any vacuum black hole characterized by its mass and angular momenta, the 
Hessian of the entropy always has a negative eigenvalue. Since the Hessian of the entropy 
is related to the inverse of the Hessian of the Gibbs potential [23] , this implies at generic 
points in the moduli space, i.e. away from the hypersurfaces defined by a vanishing eigen- 
value, an instability in the grand canonical ensemble [22J. 

Identifying which doubly spinning black hole solutions exhibit this ultra-spinning regime 
will be our next objective. We will first find the threshold of the black membrane regime for 
MP black holes with at least two large angular momenta. To compare the ultra-spinning 
regimes of both, the black hole and black ring, we present a careful analysis of their ultra- 
spinning regime. For the black ring we observe a similar but slightly different ultra-spinning 
phase. That is, the temperature does not have a minimum but rather there is a 'turning 
point', which is responsible for the boosted black string behaviour. We also find that even 
after adding the second angular momentum the ring can have a membrane-like behavior. 
However, for large enough values of the second angular momentum the 'membrane phase' 
disappears. 

The dynamical instabilities were related to the thermodynamic ones by (a conjecture 
of) Gubser and Mitra [23] : gravitational backgrounds with translationally invariant horizon 
develop a tachyonic mode (negative non-conformal mode) whenever the specific heat of the 
black brane geometry becomes negative. Our goal is not to study dynamical instabilities, 
which would require a perturbative analysis, but rather to generalize the arguments of [H] 
to black rings. However, to connect the work of [9] to [21], one has to understand the black 
hole thermodynamics in the grand canonical ensemble (see, e.g., [22]). At first sight, it is 
surprising that an analysis in the microcanonical ensemble [9] provides information about 
the membrane phase. We point out a subtle relation between the microcanonical and grand 

3 The complex geometry is obtained by the usual analytic continuation of time coordinate, r = it. 
4 For a black ring with one angular momentum the 'isothermal compressibility' (moment of inertia) is 
always negative [7]. 
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canonical ensembles, which is valid for the analytic solutions we are interested in, but it 
may be more general, and argue that this is why the microcanonical ensemble also encodes 
information about the membrane phase. 

The remainder of this paper is organized as follows. We start in Section 2 with a brief 
review of the counterterm method for asymptotically flat spacetimes. We then compute the 
stress tensor and the corresponding asymptotic charges for black objects with two angular 
momenta: Myers-Perry black hole, black ring, and black branes. In Section 3, we present 
an analysis of the thermodynamic stability of the doubly spinning ring. We compute the 
thermodynamic action and check the quantum statistical relation. We also analyze in great 
detail the response functions. In Section 4, we investigate the black string/membrane phases 
of doubly spinning black objects. Finally, we conclude with a discussion of our results. In 
Appendix A we present the expressions of angular velocities and temperature for a general 
metric with two angular momenta. Appendix B contains some general aspects of black 
holes thermodynamics, the local stability conditions, and concrete expressions for some of 
the response functions used in Section 2. 



2 Stress tensor and conserved charges 

In this section we apply the counterterm method to doubly spinning five- dimensional vac- 
uum solutions of Einstein gravity. We explicitly show how to compute the boundary stress 
tensor and the conserved charges for Myers-Perry black hole, doubly spinning black ring, 
and doubly spinning black branes. 



2.1 Quasilocal formalism and conserved charges 

To begin our considerations on thermodynamics of doubly spinning black objects in five 
dimensions, we recall the description of quasilocal formalism [12] supplemented with coun- 
terterms. 

To define the conserved charges we use the divergence-free boundary stress tensor pro- 
posed in [TT] : 

Ti3 = ~^mw^ = s^g~ 5 ( Kij ~ hijK ~ * (7 ^' _ n hij) ~ hijU ^ + (2,1) 



where \& = y/3/y/2TZ, hij is the induced boundary metric, and IZij is its Ricci scalar. A 
rigorous justification and more details about this proposal can be found in [26, 125] . 
Here, I is the renormalized action that includes counterterms, 



R^jd 5 x + -^— [ \K-J^)^\h\d i x (2.2) 



16nG 5 J M 8itG 5 j qm 

K is the extrinsic curvature of DM and e = +1(— 1) if dM is timelike (spacelike). 
The boundary metric can be written locally in ADM-like form 

h ij dx l dx j = —N 2 dt 2 + a ab (dy a + N a dt) (dy b + N b dt) (2.3) 
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where N and N a are the lapse function and the shift vector respectively and {y a } are 
the intrinsic coordinates on a (closed) hypersurface S. If the boundary geometry has an 
isometry generated by a Killing vector £ l , a conserved charge 

0£ = j> d*y^n%e (2.4) 
can be associated with the hypersurface £ (with normal ri 1 ). 

2.2 Doubly spinning solutions 
2.2.1 Black hole 

The Einstein equations in higher dimensions have spinning black hole solutions [6J. In five 
dimensions, the Myers-Perry black hole in Boyer-Lindquist type coordinates is 

ds 2 BH = -dt 2 + E dr 2 + dO 2 ^ + (r 2 + a 2 ) sin 2 9 dcf 2 + (r 2 + b 2 ) cos 2 9 d^ 2 

+ ^ (dt - a sin 2 9d<p-b cos 2 9 dip) 2 (2.5) 

where 

£ =r 2 + a 2 cos 2 # + 6 2 sin 2 #, A = (r 2 + a 2 )(r 2 + b 2 ) - mr 2 (2.6) 

and m is a parameter related to the physical mass of the black hole, while the parameters 
a and b are associated with its two independent angular momenta. This metric depends 
only on two coordinates, < r < oo and < 9 < n/2, and it is independent of time, 
— oo < t < oo, and the azimuthal angles, < (ft, ip < 2n. 

Since r is playing the role of a radial coordinate in this coordinate system, the event 
horizon is also the null surface determined by the equation g rr = . So, the event horizon 
of the black hole can be computed by using ( ]A.5j) . which implies A = 0. The largest root 
of this equation gives the radius of the black hole's outer event horizon 

r 2 = 1 (m - a 2 - b 2 + \/{m — a 2 — b 2 ) 2 — 4 a 2 6 2 ) (2.7) 

Notice that the horizon exists if and only if 

a 2 + b 2 + 2\ab\ <m (2.8) 

so that the condition m = a 2 + b 2 + 2|a6| or, equivalently, r 2 = \ab\ defines the extremal 
horizon of a five dimensional black hole (when one angular momentum vanishes, the horizon 
area goes to zero in the extremal limit). Otherwise, the metric describes a naked singularity. 
In the asymptotic limit, r^oo, the metric (12. 5p approaches Minkowski space 

ds 2 = -dt 2 + dr 2 + r 2 (d9 2 + sin 2 9 dcj) 2 + cos 2 9 dip 2 ) (2.9) 
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We use the expression of black hole metric in Boyer-Lindquist coordinates to compute the 
boundary stress tensor and we obtain the following non-vanishing components: 



3 1 5. 9 , 9 .cos20 
1 / „ sin 2 



8tt& 5 \ 2 r A 3 r 



its 



8nG 



-2am S ^- + 0(l/r 5 )) 
r 6 J 



1 (2 , , 0N cos 29 „„. , „.\ . 

(a 2 - b 2 ) — — + 0(l/r 3 ) J , (2.10) 



8ttG 5 \3 



1 12 + 

G 5 \3 v ' r 



8ttG, 

1 (2 , 2 , 2 , (1 + 2 cos 20) cos 2 ^ , „ 

1 / , cos 2 sin 2 6> B , 
7* = _^-4a6m ^ + 0(l/r 5 , 

This stress tensor is covariantly conserved with respect to the boundary metric (12. 9p . We 
also notice that, for equal angular momenta, the diagonal 'angular' components of the stress 
tensor vanish — this is intuitively expected due to the enhanced symmetry. 

Using the definition (12. 4p . it is straightforwardly to obtain the conserved charges asso- 
ciated with the surface S as 

M = jf d 3 2/vW T i:j & , = £ tfyy/an* r y f J , = £ tfyy/an* r tj f j 

where the normalized Killing vectors associated with the mass and angular momenta are 
£t = 9t 5 = ^5 an d £v> = ®ip respectively. We find 

3irm ixma 7rm& 

which is in perfect agreement with the ADM calculation. 



2.2.2 Black ring 

A black ring is a five-dimensional black hole with an event horizon of topology S 1 x S 2 
and the metric was presented in jl] - the solution of Emparan and Reall has one angular 
momentum. In five dimensions, a more general solution for a black ring with two angular 
momenta was presented by Pomeransky and Sen'kov Some properties of the solution 
including the structure of the phases in the micro canonical ensemble are discussed in [29J . 
A study of its geodesies has been performed in |30j and a careful investigation of global 
properties appeared recently in [31] . We provide here a brief account of the doubly spinning 
black ring solution and compute the boundary stress tensor and the conserved charges. 

We will use the solution in the form presented in [5]. The metric depends just on the 
coordinates x and y defined within the following intervals — 1 < x < 1 and — oo < y < — 1. 
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Notice that x is like an angular coordinate — this observation will be useful when we will 
define new coordinates that make asymptotic flatness clear. 

The metric has a coordinate singularity where g yy diverges. The event horizon of the 
doubly spinning black ring is located at the smallest absolute value of 1 + A y + vy 2 = 0, 
namely 

VH = Yv (2-12) 

For a regular black ring solution, the parameters v and A are constrained to satisfy [5]: 

< i/ < 1 , 2v^<A<l + z/ (2.13) 

In the limit v — > the black ring with one angular momentum (J^) is recovered (J^ is the 
angular momentum on S 2 ). The limit A — > 1\fv was carefully studied in (29] and shown to 
correspond to regular extremal black rings. 

We use a coordinate transformation similar to the one in 



x = -l+4k 2 a 2 — y = -l-4k 2 a 2 — =— , a=\ 2.14) 

r z r z V 1 — v 



In these coordinates d t , d^, and are Killing vectors and the asymptotic metric is the 
same as (12. 9ft . 

The boundary stress tensor in these new coordinates is 

1 / 12fc 2 A 1 8g\M cos 20 5 , 



8ttG 5 V. (l + ^-A)r 3 3(l + i/- A)(l - z/) 2 r 

1 / 16fc 3 A (1 + A - 6v + \v + z/ 2 ) ^/(l + z/) 2 - A 2 sin 
T * = 8^1 (l + z,-A) 2 (l-z,) 2 ^ + 



1 / 32^^(1 + z,) 2 - A 2 ] cos 2£ _ 

r n = —r —p-t — — + O 1/r- 

87rG 5 I (1 + z/ — A)(l — v) z r d 

1 / 2k 2 F 2 [u,X] cos 20 , ., 



0(iy)j , (2.15) 
+ C(l/r 3 



8ttG 5 V3(l + ^-A)(l-z/) 2 r 

1 fk 2 (-F 3 [u, A] + F 4 [u, A] cos 20) sin 2 
8irG 5 V 3(l + z/-A)(l-z/) 2 ~ 
1 / fc 2 ^M+%j cos 20) cos 2 _ 3 
^ ~ 8ttG 5 V 3(l + v-\)(l-v) 2 r U[/r 

1 / 128fc 4 Ay^(4Az/ - (A 2 + (1 - z/) 2 )(l + z/)) cos 2 sin 2 5 



where 



Fi [z/, A] = 1 - 5z/ - z/ 2 + 5z/ 3 + A 2 (3 + 7z/) + A(l - 14z/ - lv 2 ) , 

F 2 [u,X] = l-llz/-z/ 2 + llz/ 3 + A 2 (3 + 13z/)+4A(l-5z/-4z/ 2 ), 

F 3 [v,\] = 5-7z/-5z/ 2 + 7z/ 3 + A 2 (15 + 17z/) -4A(l + 13z/ + 2z/ 2 ), 

F 4 [z/,A] = 7-29z/-7z/ 2 + 29z/ 3 + A 2 (21 + 43z/) + A(4-92z/-40z/ 2 ) 
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As in the case of doubly spinning black hole, this stress tensor is covariantly conserved 
with respect to the boundary metric ( 12. 9ft . However, since for the doubly spinning black 
ring the angular momenta can not be equal, namely 




(2.16) 



there is no similar symmetry enhancement as in the black hole case in the angular part. 

By plugging the expressions of the boundary stress-energy components (12. 15ft in ( 12. 4j) 
we find the following expressions for the conserved charges: 



As expected, the charges computed by using the quasilocal formalism recover correctly 
the ADM results [5]. 

In principle, one can obtain a black hole and a black ring with the same conserved 
charges. However, an asymptotic observer can not distinguish between a black hole and a 
black ring just by computing the conserved asymptotic charges. We would like to emphasize 
that it is expected that the subleading terms of the quasilocal stress tensor encode the 
information necessary to distinguish between black objects with different horizon topologies. 

2.2.3 Black membrane 

Here we would like to apply the quasilocal formalism to doubly spinning black p-branes, 
dubbed also black membranes (BM) or black strings if p = 1. The black membrane metric 
we are interested in is obtained by adding flat directions to a 5-dimensional black hole with 
two angular momenta. Therefore, the metric is 



where ds 2 BH is the black hole metric defined in (I2.5p . 

Since the number of dimensions and the topology are changed, one expects changes with 
respect to the former discussion. For example, the form of the counterterm leading to a 
finite actions may be different when the number of dimensions is increased. However, in 
this particular case, what is important is the 'seed' 5-dimensional solution to which we add 
the flat directions. Thus, the form of the counterterm does not change but the stress tensor 
will have new components. 

A similar computation as for the doubly spinning black hole reveals that the stress tensor 
of the BM is the one in (I2.10p supplemented with the components in the new directions: 





(2.19) 




(2.20) 



This result resembles the tension (per unit length) of the black string. 
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3 Thermodynamic instability of the black ring 



In this section, we discuss the thermodynamics of a doubly spinning ring in the grand 
canonical ensemble. 

So far, we have computed the conserved charges of neutral spinning black objects with 
two angular momenta by using the quasilocal formalism. However, the quasilocal formalism 
is a very powerful tool for understanding the thermodynamics in more detail. In particular, 
one can compute the action and, therefore, the thermodynamic potential. 

In what follows, we present a detailed analysis of thermodynamic stability of the doubly 
spinning black ring — an analysis of the thermodynamic stability of Myers- Perry black hole 
with two angular momenta can be found in |2U] . 

Let us start by computing the angular velocities and the temperature for this solution. 
From (IA.6j) we obtain the following expressions for the angular velocities: 



A (l + „) _ (1 _ jjygTg 1 + U-X 1 1 + U-X f , 

Vi + ^ + A' * ~ 2k\TTVTx {3 ' h 

The area of the event horizon and the temperature (1A.7|) are 



_ 32vr 2 A; 3 A(l + A + z/) yg^4gl - u)(y h 1 - y h ) 

(l- u )2(y-i- Vh ) ' 8vrA;A(l + A + ^) 1 " ] 

Note that yh = ~ A+ ^ 2 ~ — is the biggest root of (1A.5jl which corresponds to the outer event 
horizon — at this point, it might be useful to emphasize again that — oo < y < — 1. 

The starting point of the Euclidean approach to black hole thermodynamics is the 
partition function |16^l 

Z(f3)= [ d[<7,0]e-« (3.3) 



where is a collective notation for the matter fields, d\g, <fi] is the measure, and I\g,<j>) is 
the Euclidean classical action. The gravitational partition function is defined by a sum over 
all smooth geometries (including black holes) that are periodic with period /3 = T -1 in the 
same class of boundary conditions e.g., AF spacetimes. 

For our purpose it is enough to consider the saddle point approximation. The grand 
canonical partition function is then Z = Tre~ /3( - H ~ naJa ^ ~ e~ Icl (here we are interested in 
black objects with two angular momenta), where I c i is the classical action. The saddle 
point is usually referred to as a gravitational instantonu 

The thermodynamic (effective) potential associated to grand canonical ensemble is 

G[T, Q a ] = k = M-TS-Q a J a (3.4) 



5 It should be understood as a low energy effective theory rather than a proper theory of quantum gravity. 

6 A quantum field can be treated as a small perturbation about the gravitational instanton. The next 
order contribution, which gives the one loop correction, includes also the thermal radiation outside the 
black hole. 
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On the Euclidean section, the topology near the horizon is modified^ and one has to 
deal with manifolds with conical singularities. It was shown in (32J S3] that the conical 
defect has a contribution to the curvature and, consequently, the path integral is rescaled 
by e s . However, this can be intuitively interpreted as a consequence of a trace over the 
macroscopically indistinguishable microstates. 

Let us now compute the action for the doubly spinning black ring. Since the Ricci scalar 
vanishes on-shell, the only contribution to the action is coming from the surface terms. To 
evaluate these terms, it is convenient to use the (r, 6) coordinate system instead of the (x, y) 
coordinates — the reason is that the normal to the boundary has just one non- vanishing 
component. We find 

\ = 2 ^(A (1 -^-^,Ai ro ^ )5i , l2f) + o(i/r) (35) 

where F 5 [u, A] = 1 + 3 A 2 + 6u + 5u 2 — 4A — 8Az/. The expression for the total action is 

nk 2 A 

and satisfies H3.4j) . which is the quantum statistical relation for the doubly spinning black 
ring. This can also be regarded as a non-trivial check that the entropy S of this solution 
is, indeed, one quarter of the event horizon area Ah- 

We have checked that the usual thermodynamic relations 

are satisfied and so the Gibbs potential G[T, f^, fl^} is indeed the Legendre transform of 
the energy M[S, J<p, J^} with respect to S, J^, and J^. 

We want to also point out that, in the light of the new developments in understanding 
the balance condition for gravity solutions [T5], the form of quantum statistical relation 
hints to the fact that this solution is balanced. Indeed, our results are in perfect agreement 
with the recent detailed analysis of the global properties of the doubly spinning black ring 

Now, we are ready to discuss the thermodynamic stability in the grand canonical ensem- 
ble — in Appendix [B] we summarize the thermal stability conditions and present explicit 
expressions for some response functions we are interested in. We analyze in detail the 
response functions that signal the (in) stability of the black ring against fluctuations. 

We consider first the specific heat at constant angular velocities 



Cn^T(^) (3.8) 



\dT 



7 The origin in the Euclidean spacetime translates to the horizon surface in the Lorentzian spacetime. 
The Euclidean section can be understood as an effective description where the microstates can not be 
distinguished. 
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Figure 1: Scatter plots in parameter phase space (u, A) for the doubly spinning black ring. The 
plot on the left shows the regions (10,000 points) where the heat capacities are negative, Cq < 
(gray) and Cj < (black). The regions where the compressibility (gray) and the det[e] (black) 
are negative cover the entire parameter space (plot on the right) implying the local thermal 
instability of the doubly spinning black ring. The region in the parameter space is bounded: 
< v < 1 and A by the functions 1 + v and shown as the dashed and solid lines respectively. 

The analytic form of this quantity is too complicated to be written down here. Instead, we 
show on the left hand side in Fig. [T] a scatter region in the parameter space of the doubly 
spinning black ring where this heat capacity Cq is negative (gray - 10,000 points). Note 
also that the parameters in the solution (12.131) are constrained and represented as a dashed 
line for A = 1 + v and solid line for the extremal black ring with A = 2-s/u. 

In a similar way, we explore the region where the specific heat at constant angular 
momentum 



is negative. The region (in black) where Cj < is shown in the scatter plot on the left of 



We observe a region in the parameter space of the doubly spinning black hole where both 
specific heats, Cq and Cj, can be positive simultaneously. However, this condition is not 
sufficient to draw the conclusion of thermodynamic stability: one should also investigate 
the matrix of 'isothermal moment of inertia'. 

These response functions are defined as 



We observe in Fig. [T] that the spectrum of the matrix of isothermal moment of inertia, 
spec[e a fe], is nowhere positive definite in the parameter space. 

Since there is no overlap region in parameter space in which all the response functions 
of interest are positive definite, we conclude that the doubly spinning black ring is unstable 
in the grand canonical ensemble. 




(3.9) 



Fig. m 




(3.10) 
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4 Instabilities from thermodynamics 



Many known stationary black holes in higher dimensions present a black string or, more 
general, black brane phase — we will refer to it as the 'membrane phase'. That is, as the 
angular momenta are sufficiently increased (the ultra-spinning regime), the behaviour of 
some black holes and black rings changes to that of extended black branes and strings. 

In the next subsection , we deal with the ultra-spinning black holes. From the study of 
the Gibbs potential's Hessian, we show the existence and find the locus of the transition 
points to the membrane phase. We also argue that there is a subtle relation between the 
microcanonical and grand canonical ensembles that may be at the basis of some of the 
results for ultra-spinning black holes discussed recently in [TP] . 

The analysis can be extended to (doubly) spinning black rings [§. These results are 
presented in Section 14.21 

To compare different examples we will make use of the dimensionless expressions for 
the temperature t, the spin j, and the angular velocity u defined by 

t D ~ 3 = c t GMT D ~ 3 u D -* = c^GMn D -\ i D ~ 3 = c i-Gj^2 
where the numerical constants are 

2 (4vr)^ /D-3\^ 16 p-3)^ _ Cg ~ 2) D - 2 U9 , 

Ct (D-2) D _ 3 \D-A) ' Cu (D-2) , 3 2W (fl _ 3) V' 1 J 



4.1 Ultra-spinning black holes 

Due to the qualitative changing behaviour of black holes as the dimensions are increased, the 
authors of [9] have argued that the ultra-spinning black holes — those in D > 6 dimensions 
that can have arbitrary large angular momentum per unit mass [6] - become unstable. 
The transition of these black holes, from behaving like a spherical black hole to behaving 
like a black membrane as the spin grows, was established to be at the minimum of the 
temperature. From that point onwards, the temperature increases in a similar way as for 
the black brane temperature. 

The minimum of the temperature where the behavior of the singly spinning black hole 
changes is determined as [9] 

a 4 = (4-3) 
r z h D — 5 

This result was also obtained by using a different method, namely finding the divergences 
of the 'Ruppeiner curvature' |23 | . It was shown in [35] that, for a singly spinning Myers- 
Perry black hole, this curvature^ blows-up exactly at the value (14. 3[) signaling a thermal 
instability of the system. 



8 The doubly spinning solutions we consider are with the spins in orthogonal planes. Other black hole 
solutions, where the two spins are parallely oriented, were also studied [34] . 

9 Thc Ruppeiner curvature is the scalar curvature of the Hessian matrix of the entropy. 
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A qualitative understanding of this fact is related to the observation that, as the spin 
becomes large, the event horizon spreads out in the plane of rotation: it becomes a higher 
dimensional 'pancake' approaching the geometry of a black brane. 

The existence of the ultra-spinning limit resembling black branes has a remarkable 
consequence. Black branes were shown to be classically unstable [11] so that the ultra- 
spinning black holes would inherit the Gregory-Laflamme instability. The threshold of the 
classical instabilities and the connection to the thermal instability as conjectured by [24J 
(see, also, [25]) requires a linearized analysis of the perturbations of the black hole solutions. 

However, the transition to a membrane-like phase of the rapidly spinning black holes 
can be established from the study of the thermodynamics of the system. The existence and 
location of the threshold of this regime is signaled by the minimum of the temperature and 
the maximum angular velocity as functions of the angular momentum. 

It was observed in [10] that, for ultra-spinning black holes, this is in tight correspondence 
with a vanishing eigenvalue of the Hessian of the Gibbs potential. A complete thermody- 
namic analysis, though, should be based on the full Hessian of the thermodynamic potential 
rather than only a s tndy of the determinants We will see in the next subsection that the 
membrane phase of a doubly spinning ring is not signaled by a zero-eigenvalue of the Gibbs 
potential's Hessian. 

For ultra-spinning black holes, there is a direct relation between (some response functions 
in) the microcanonical and grand canonical ensembles§ To see that, let us compare the 
expressions of two particular response functions in these two ensembles: 

fd 2 s\ i fdsi\ si (dr\ , fS 2 G\ (aj\ 

W) M = -T\jj) M + ^{ej) M and {«e) T = -{aa) T (4 ' 4) 

We have checked that in the particular case of the singly spinning black hole, indeed, 
these two response functions are inverse proportional at the particular point where the 
temperature has a minimum. Therefore, an inflexion point in the microcanonical ensemble 
corresponds to a divergence of the corresponding response function in the grand canonical 
ensemble. This may well be an explanation for the results obtained in [10]. Moreover, this 
point should not be considered as a sign for an instability or a new branch but a transition 
to an infinitesimally nearby solution along the same family of solutions. The numerical 
evidence of [lU\ supports this connection with the zero-mode perturbation of the solution. 

We now examine the situation for a more general family of ultra-spinning Myers-Perry 
black holes with multiple spin parameters, a iy where % = 1,2,...,N and N = [(D — l)/2]. 
The black hole is characterized by the mass parameter [i and the horizon radius (the 
largest root of) 

1 N 

A* = T+5 II(^ + °*) ' ( 4 - 5 ) 

r h i= l 

A spinodal is defined as a line separating the regions of stability and instability of a homogeneous 
system. It is important to emphasize that all spinodals are zero-determinant lines, but in general not all 
zero-determinant lines are spinodal. 

11 Different ensembles correspond to different physical conditions and so, in more general cases, one does 
not expect such a relation. 
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by which we can express the thermodynamics 

M = ifi n D ~ 2 ^' ~\ a n ^ = 2 , 2 » 

^ = ^- 2 /ir,, T=^- (r^-j-L_ , (4.6) 



=i 



where e = moc^-D- A sufficient, but not necessary, condition for the existence of ultra- 
spinning black holes was given in [§]. In even(odd) dimensions at least one(two) of the 
spins should be much smaller than the rest. The ultra-spinning regime is obtained in the 
limit 

< ai, a 2 , a k « a k+ i, a N ->■ oo (4.7) 

where N — 1 > A; > 1 + e. The generic limiting black brane metric whether static, with all 
finite angular momenta ai,..,a k vanishing, or spinning, with some a\,..,a k non- vanishing, 
is the product 5-^-2 (;v-*+i) x M 2(v-fc)_ 

Our focus will be on the case in which the black hole has at least two large spins and 
we set the remaining angular momenta to zero. When the angular momenta are equal, 
Jk+i = ■■■ = Jn = J , the Ruppeiner curvature scalar blows up ato 

a 2 D-3 



i 2/c-l 



According with the arguments in [23J, this signals a thermodynamic instability. However, 
the expected new phase should correspond to the black membrane phase of ultra-spinning 
black holes and not to a new branch of solutions. 

This is further supported by examining the eigenvalues of the Hessian of the Gibbs 
potential. Indeed, we find that the divergences of the Ruppeiner curvature pinpoint the 
zero of the determinant of the Gibbs potential's Hessian. 

Also, by studying the temperature 



(D - 3) (l 



4J 2 



(D-3) S 2 



AS^ (1 + ^)^ 



D-l + n 

2) 



n = 2k-l-e (4.9) 



we find that the temperature has a minimum at exactly (14.81) . while the angular velocity 
Q reaches its maximum value. Therefore for these more general ultra-spinning black holes, 
similarly to the singly spinning situation discussed in [9], once the minimum is reached 
the temperature increases and the angular velocity decreases signaling a transition to a 
membrane phase. This conduct is shown in figure [2] (J points on the solid thin line) for the 
particular case of D = 7, k = 2 so j\ = j'2 = 3<j> and j'3 = 0. 

Another case of interest is the ultra-spinning black holes that resemble spinning black 
branes, when some of the slower spins are non- zero a±, a k ^ 0. It is not our goal to make 
a detailed analysis of this case here. Nevertheless, in all the cases where the non- vanishing 
spins are set equal, we find divergences of the Ruppeiner scalar curvature which could help 
to detect the threshold of their membrane phase. 



12 Note that ^ = H^- and in the particular case when k = k max our results agree with those of [35] . 
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Figure 2: Plots of the temperature (left hand side) and angular velocity (right hand side) as 
functions of the angular momentum j^, for a fixed mass, for different black objects. These include 
the singly spinning Myers-Perry black holes in five dimensions of space-time (black dashed line) and 
its seven dimensional cousin (solid thin line). The singly (solid thick) and doubly spinning black 
ring (light gray) for different values of angular momenta (right towards left) = 0.1,0.2 on the 
S 2 , are also shown here. For this same values of the second angular momenta the five-dimensional 
doubly spinning black holes are represented by the dashed light gray lines. 



4.2 Membrane phase of black rings 

Other solutions, e.g. the black ring with one angular momentum, also exhibit an ultra- 
spinning behavior. The black ring, which is characterized by the radii r and R of the 
spheres S D ~ 3 and S 1 , respectively, becomes thin in this limit (when r << R)E 

Since the final expressions for the response functions are very complicated for the doubly 
spinning ring (see Appendix B), we prefer to present the 'conjugacy diagram' of the angular 
velocity versus the angular momentum and the plot of the temperature as a function of the 
angular momentum, both for a fixed mass (see Fig. [2]). 

For the singly spinning black ring, an analysis of the temperature as a function of the 
angular momentum was presented in |38j. In this case (solid thick line in the plot on the 
left hand side of Fig. 2), the temperature does not have a minimum, but there exists a 
turning point In our analysis, the turning point II for the black ring plays a similar 
role as the minimum of the temperature for the black hole. That is, it signals a change 
in the thermodynamic behaviour of the black ring. In fact, it is the starting point of the 
ultra-spinning regime where the black ring can be approximated by a boosted black string. 

Using the Poincare 'turning point' method, this special point was carefully studied in 
[38J. In particular, they found a divergence of the Ruppeiner curvature. In the conjugacy 
diagram (on the right) there is also a turning point II at the same minimum value of the 
angular momentum j^. 

The question is then if there still is a relation between the microcanonical and grand 

13 This thin regime was essential to find perturbatively the higher dimensional black rings cousins. More- 
over, a generalization of this construction to black branes led to the construction of blackfolds [551 ISZ] • 
14 We would like to point out that in (Sherk-Schwarz-)Anti-de Sitter, there is also a turning point [39] . 
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canonical ensembles in this case. We have explicitly checked, using the results of [7J, that 
one of the eigenvalues of the Hessian of the Gibbs potential is zero at this specific point II 
while the second eigenvalue never changes its sign. Therefore, we conclude that the turning 
point is the onset of the ultra-spinning black string phase. 

A far more richer structure is found for the doubly spinning black ring. The angular 
momentum on S 2 is bounded as G [0, 1/4] and for a specific the black ring can always 
be extremal (in the limit A — > l^fv as shown in [29] ). But besides extremality, according 
to how large is, the behaviour of the doubly spinning black ring changes. There are 
two distinctive regions in the microcanonical ensemble. On one hand, for < < 1/5, 
there are phases with the characteristic cusp for black rings with the two (fat and thin) 
branches. But on the other hand, for 1/5 < < 1/4, the fat black ring branch disappears 
and so there are no cusps. As we will discuss in what follows, this will become relevant to 
understand the physics and regimes of the doubly spinning black ring. 

To explore how the physics of the black rings at fixed mass is modified as we turn on 
the angular momenta along the S 2 , j^, we will study the temperature and angular velocity 
as functions of the S 1 angular momentum, j^, for different fixed values of j^. 

Doubly spinning black rings with < < 1/5 

In this case the situation is similar, to some extent, to what we found before for the singly 
spinning black ring. In this range, as we turn on j^, there also are turning points with 
tangents of infinite slope signaling the onset of the black membrane phase that coincide 
with the cusps (in parameter space), namely at A = — (1/4) (1 + v — (9 + z/) 1 / 2 (l + 9//) 1 / 2 ). 
Figure [2] shows this change in behavior explicitly. On the left (light gray curve) in the t vs. 
jtl, diagram, we observe the turning point 77/ that corresponds to the minimum value of 
jtj) angular momentum. In the conjugacy diagram uj vs. (also the light gray curve), the 
corresponding point (III) is a turning point. 

But an interesting difference with the single spin black ring occurs for a large enough 
S 1 angular momentum. The temperature of the black ring in its membrane phase increases 
while the area of the event horizon decreases up to a point where the spin-spin interaction 
is large enough making a turn to abruptly become extremal with zero temperature. This is 
the maximum critical value labeled IV in figure [2j Therefore, the black membrane phase 
exists between points III and IV. 

Finally, as for the singly spinning black ring, we computed the eigenvalues of the Hess[G] 
and found that none of them are zero the relevant turning points. 

Doubly spinning black rings with 1/5 < < 1/4 

This black rings (lightest gray line in figure [2]) with larger S 2 angular momentum show no 
turning points and therefore have no membrane phase. This mimics its behavior in the 
phase diagram of microcanonical ensemble where the cusp and fat branch of these doubly 
spinning black ring disappear [29J. The lack of a fat black ring branch seems to coincide 
with the lack of a black membrane phase. Therefore such solutions would never be captured 
with long distance effective approaches [36| [37]. 
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Note that for certain fixed € [1/5, 1/4] the temperature grows, reaching a maximum 
at V and rapidly decreasing to zero to become extremal. It would be interesting to explore 
the physical meaning of these points which we observe to correspond to an inflection point 
d 2 S/d 2 J = for fixed mass. 

In summary, we have found examples for which the zero eigenvalues of the Hessian of 
the Gibbs potential can also be turning points (with tangents of infinite slope) and not 
just critical points as for the Myers-Perry black holes. Other less symmetric solutions, 
such as the doubly spinning black ring, do not show this connection between the zeros 
of the Hess[G] and the onset of the black membrane phase. Moreover, we showed that 
certain doubly spinning black ring, those with 1/5 < < 1/4, have no membrane phase. 
Therefore, these particular solutions fall into the same category as other black holes with no 
membrane phase as the four dimensional Kerr black hole and the five dimensional Myers- 
Perry black hole. 

5 Discussion 

In this paper we have analyzed in detail the thermodynamic stability of neutral doubly 
spinning black objects. We have analytically computed the response functions and pre- 
sented strong evidence showing that the doubly spinning black ring is thermodynamically 
unstable. That is, there is no region in the parameter space in which all the response func- 
tions are positive definite. 

We have provided an explanation of why the microcanonical and grand canonical en- 
sembles for ultra-spinning black holes are related in a very specific way. An inflexion point 
in the microcanonical ensemble corresponds to a divergence of the corresponding response 
function in the grand canonical ensemble. We will comment on the validity of this argument 
and on the significance of these results in the last part of this section. 

The onset of a membrane phase of different doubly spinning black objects was identified. 
We have found that the onset of the black membrane phase for all black holes that we have 
studied (except for the less symmetric doubly spinning black ring) is characterized by at 
least one zero eigenvalue of the Hessian of the Gibbs potential. A tight relation with the 
classical perturbations, where a transition to an infinitesimally nearby solution of the same 
family branch happens, is expected in the cases where the connection between membrane 
and zero-eigenvalue of HessfG] exists. The numerical evidence of [TQl [22] supports this 
connection precisely with the zero-mode perturbation of the solution. 

We now discuss the thermodynamics of doubly spinning black rings in the grand canon- 
ical ensemble. An analysis of this solution in the microcanonical ensemble was presented 
in [29]. In general, for black holes, the entropy is used to obtain the phase diagrams in 
the microcanonical ensemble while the mass is kept fixed. However, in general relativity it 
makes more sense to use the total energy instead of the entropy. The reason is that this 
would require appropriate boundary conditions. 

It is important to emphasize that it is not known how the background subtraction 



17 




Figure 3: Plot of the response function e^, as a function of A. The gray curve corresponds to the 
compressibility of the singly spinning black ring, namely v = 0. As the angular momentum along 
the S 2 is increased (dashed line left towards right) the isothermal moment of inertia for different 
values of v = 0.15,0.3,0.42,0.48,0.53 changes and becomes positive for values of v > 0.46690042 
and A > 1.40685 (shaded gray region). 



method can be applied to black rings, because it is not clear how to choose the background 
solution. We therefore used the counterterm method to compute the action and the grand 
thermodynamic potential (which is a Legendre transform of the energy). Since the expres- 
sions of the response functions are too complicated for analytical treatment (see Appendix 
B), we have plotted the regions in the parameter space where they are positive. 

Let us know compare some of our results with other well understood examples — at 
this point we are interested just in the thermodynamic instabilities, not in the relation with 
dynamical ones. The fact that the Schwarzschild black hole has a negative heat capacity 
means that the thermodynamic ensemble is dominated by diffuse radiation states rather 
than black holes states but is classically dynamically stable at the linearized level. When 
adding angular momentum the situation changes and the heat capacity can be positive 
for a large enough angular momentum. However, this condition is not enough to conclude 
that the system is thermodynamically stable. The stability also implies that when angular 
momentum is added to the system the angular velocity goes up. 

For a black ring with one angular momentum the heat capacity can be positive, but 
the momentum of inertia is always negative [7]. Therefore, the singly spinning black ring is 
thermodynamically unstable. As in the case of one angular momentum, the heat capacity 
of a doubly spinning black ring can be positive in some region of the parameter space. 
However, there is a key difference when the second angular momentum is turned on. That 
is, the component of the momentum of inertia associated to S 1 of the black ring can become 
positive — this is explicitly shown in Figure [31 

Since there are two angular momenta one should also investigate the effect of coupled 
'angular' inhomogeneities. A careful study of the determinant of the momentum of inertia 
matrix shows that there is no region in the parameter space with the desired properties and 
so the doubly spinning black ring is also thermodynamically unstable. 
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Figure 4: The microcanonical phase diagram - entropy, s, as a function of the angular momentum, 
j - for fixed mass (on the left) of the singly spinning Myers-Perry black hole in D = 5, 6, 10 
dimensions of space time ( black, gray and lightgray lines respectively) . The grand canonical phase 
diagram (on the right) - Gibbs potential, g, as a function of the angular velocity, oj - for fixed 
temperature of the black hole. The points A,B correspond to a/rh = \/3, y| for six and ten 

dimensions, and in general for D > 5 the points are a/r^ = (D — 3)/(D — 5); the change in 
convexity of the entropy corresponds to a blowing up of the convexity of the Gibbs potential. The 
point C instead, at a/r^ = and where the Gibbs potential has a change in convexity, has no 

analog in the microcanonical scheme. 

We would like now to discuss in more detail some of the results for ultra-spinning black 
holes presented in Section 4. In Fig. HJ we show the points (A,B) in the grand canonical 
ensemble that correspond to inflexion points (A,B) in the microcanonical ensemble. This 
can be quantitatively understood by comparing the particular response functions in eq. 
(14. 4 p at a very special point in the parameter space, namely where the temperature has a 
minimum. 

We have explained that this argument applies to this particular case but not in general. 
A counterexample is the 5-dimensional black hole with one angular momentum. In this 
case there is no relation between ensembles in the sense that there is no special point in 
the microcanonical ensemble which corresponds to the inflexion point (C) in the grand 
canonical ensemble. Moreover we have checked and there are no points where an eigenvalue 
of the Hessian of the Gibbs potential vanishes. Therefore, it should not be considered as a 
sign for a membrane phase — most probably, it is similar with the Schwarschild black hole 
example for which the thermodynamical instability is not related to a dynamical one. 

One can also consider the ultra-spinning black holes with some of the finite angular 
momenta non-zero. In odd(even) spacetime dimensions, the metric of an ultra-spinning 
black hole with all but two (one) of the spins finite and non-zero will reduce to that of a 
spinning black brane. As we have shown in Section 2, the counterterm method can also 
be applied to spinning black branes and the results are similar to the ones for the 'seed' 
spinning black hole solution. We have computed the renormalized action to find the Gibbs 
potential and we expect similar thermal instabilities as for the corresponding black holes. 
As we already emphasized, though, in all these examples we expect that the thermodynamic 
instabilities do not signal a dynamical instability or a new branch, but rather a transition 
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to an infinitesimally nearby solution along the same family of solutions [10J. 

It is remarkable that our study of thermodynamic instabilities provides to some extent 
information about the zero- mode of the 'Gregory-Laflamme instability' and it may well 
be the starting point for a more detailed study of dynamical instabilities. The extension 
of the dynamical stability studies to spinning black branes has not been yet developed. 
The analytic theory of perturbations is much more involved. However, we expect that the 
spinning black branes suffer from similar instabilities as the static ones. We hope that the 
observations made in this paper will be useful in future investigations of the perturbations 
of higher- dimensional spinning black rings, black holes and spinning black branes. 
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A Temperature and angular velocities 

Consider a general stationary 5-dimensional metric that corresponds to a black object with 
two angular momenta^ : 

ds 2 = g t t(x)dt 2 + 2g t( j ) {x)dtd(j) + 2g t ^(x)dtdtp + g$${x) d(p 2 

+ 2g^p(x) d<pdijj + g^(x) dtjj 2 + g a p{x) dx a dx 13 (A.l) 

d t , d^, and are Killing vectors. Rewrite the metric in the ADM form 

ds 2 = -N 2 dt 2 + 7ii (dx i + N* dt) (dx j + N j dt) (A.2) 

with lapse function 

N 2 = -g tt + 9h ( N *) 2 + 9i>i> (^) 2 + 2 944 N+N* (A. 3) 

and shift vector 

jycj, _ 9t>p9H} ~ 9ipip9tcj> ^ _ 9t<f>9<j,ip ~ 9<j>4>9t-i> ^ ^ 
9%p - 9<M>9iH> 9% - 9<t><t>9i>i> 

The event horizon is obtained for 

N 2 = (A.5) 

In other words, it is a Killing horizon of d t + Sl^d^ + Q^d^, where and are the 
angular velocities defined as the shift vectors at the horizon: 

^ = -^| ff , n^=-N% (A.6) 



15 A similar analysis for one angular momentum can be found in 1 
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Black holes are thermodynamic objects: the causal structure of spacetime can influence the 
physics of a quantum field. The vacuum fluctuations near the event horizon cause the black 
hole to emit particles with a thermal spectrum. The Euclidean regularity at the horizon is 
equivalent to the condition that the black hole is in thermodynamical equilibrium. 

By a straightforward computation one can eliminate the conical singularity in the Eu- 
clidean section, (it,r), to obtain the periodicity of the Euclidean time. In this way, one 
obtains the following expression for the temperature of the black hole: 

T (N 2 y 

4ir^/g rr N 2 

We have used these definitions to compute the corresponding physical quantities of the 
doubly spinning ring. 



(A.7) 

H 



B Conditions for thermodynamic stability 

In this appendix, we present the conditions for the thermodynamic stability and we also 
give some useful explicit expressions for the response functions used in Section 4 — we 
follow closely |41j . 

For simplicity, let us start with a black hole with one angular momentum. We are 
interested in the thermodynamic potentials: the energy and its Legendre transforms. 

The basic extremum principle of thermodynamics (for the entropy S) implies both that 
dS = and that d 2 S < 0. The second condition determines the stability of predicted equi- 
librium states. The stability criterion in energy representation requires that an equilibrium 
state at fixed S and J is a state of minimum energy, namely a minimum of E[S, J]. The 
local stability conditions ensure that inhomogeneities of either S and J separately 

d 2 E\ fdT\ fd 2 E\ fdQ\ , x 

^ '-Z5 >0, _ = _ > (B.l) 



OS 2 J j \dSJj- ' \dJ 2 J s \dJ J s 
and also that a coupled inhomogeneity of S and J together 

f)2 rp o2 771 / 02 771 \ 2 

det(Hes 8(B ))^^-(— ) < B ' 2 > 

do not decrease the energy. 

In more generality the stability criterion states that the thermodynamic potentials are 
convex functions of their extensive variables and concave functions of their intensive vari- 
ables (see, e.g., [4"I]). 

For a grand canonical ensemble defined at fixed temperature T and angular velocities fl a 
(intensive variables) the associated potential, the Gibbs free energy, satisfies the following 
relations: 

G[T, tt]= E -TS - 07 , dG = -SdT - JdVL (B.3) 
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In this case, the local stability conditions following from the convexity of the Gibbs function 
yield 

d 2 G\ fdS\ fd 2 G\ fdJ\ .„ x 

" »F <°> «FvT =" «o <° ( B - 4 ) 



dT 2 y n \dTJ n - ' \dtt 2 J T \dnj 



T 



and 



«^^^y >0 (B.5) 

Equivalently, the heat capacities (Cq, Cj) and the 'isothermal moment of inertia' or 'com- 
pressibility' e = (dJ/dVt) T should be positive definite. 

A generalization for two angular momenta is straightforward (see, also, |20j). The 
Hessian is a 3 x 3 matrix 



Hess(G) = (-1) 



tab 



where the matrix components are a a = (^) n , tab = ffnM 5 an d the indices cover the 
angular directions a,b = <f>, if). 

Considering the relationship between the specific heats Cn = Cj + T (e -1 ) a & a a a b it 
can be shown that a thermodynamically stable system is characterized by positive heat 
capacities Cq > and Cj > and, also, a positive definite matrix of isothermal momenta 
of inertia, i.e. spec[e a ;,] > 0. 

The ((/>0)-component of the isothermal moment of inertia tensor is 

4fc 4 7rA(l + A + v) F(X,u) 

eH ~G 5 (is - 1) 4 (1 - A + z/) 2 TA^^MA 2 + X 2 v - 8u) + A 3 (l + u) - 2X(Au + v 2 - 1) 
where 

F(\,v) = A 6 (l + vf + A 5 (l + v) 2 (\ + VA 2 - Tv + v) 

+A 4 (1 + 1/) [4 + VX 2 - Av + i/(2(l + VX 2 - Av) + i/(-24 + VA 2 - 4i/ + 2i/))] 

-16VA 2 -4i/ + + 14 + i/(18 + (-14 + z/)z/))] 

+2A 3 (1 + 1/) [2 + VA 2 - 4z/ + 25 + VA 2 - 4z/ + i/(13 - llVA 2 -4z/ 



+z/(-7 + VA 2 - 4z/ + 1/)))] + 2A 2 (1 + i/)[2 + VA 2 - 4i/ + 32 - 26VA 2 - 4z/ 
+z/(32 + 14V A 2 - Av + i/(16 - 6VA 2 - 4z/ + (-2 + VA 2 - Av)v)))} 
-4A[-1 + i/(10 + 4VA 2 - 4z/ + 89 - 12VA 2 - Av + i/(-12(-6 + VA 2 - 4z/) 
+z/(-23 + 4VA 2 - 4z/ + (-2 + z»)))] 
The determinant is 



4fc 8 7r 2 A 2 (A - VA 2 - Av){X + VX 2 - Av) 2 ^{l + A + vf G{X,u) 
(G 5 ) 2 (-l + z/) 4 [-A 2 + (1 + u) 2 f/ 2 [u(-X 2 + (1 + v) 2 )f/ 2 Z(A, 1/) 



(B.6) 



where 



Z(X,v) = 8VA 2 - Aw - A 3 (l + v) - AVA 2 - 4i/(l + v) + 2A(-1 + 4i/ 



.2\ 



;/ 



z/ 2 l 



G(A,i/) = A 5 + 7A 4 (1 + 1/) - A 2 (l + + 3VA 2 - 4z/ + (2 - 3VA 2 - 4i/)i/ 

+A 3 [5 - 6VA 2 - 4z/ + (26 + 6VA 2 - 4z/)z/ + 5z/ 2 ] - 8i/(l + Hi/ + llz/ 2 + z^ 3 ) 



-A[(8 - 27VA 2 - 4z/)z/ + 9(16 + 3VA 2 - Av)v 2 + (8 + 3\/A 2 - 4z/)z/ 3 - 3 a/A 2 - 4z/] 
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